Driven translocation of a polymer: fluctuations at work 
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The impact of thermal fluctuations on the translocation dynamics of a polymer chain driven 
through a narrow pore has been investigated theoretically and by means of extensive Molecular- 
Dynamics (MD) simulation. The theoretical consideration is based on the so-called velocity Langevin 
(V-Langevin) equation which determines the progress of the translocation in terms of the number 
of polymer segments, s(t), that have passed through the pore at time t due to a driving force /. 
The formalism is based only on the assumption that, due to thermal fluctuations, the translocation 
velocity v = s(t) is a Gaussian random process as suggested by our MD data. With this in mind 
we have derived the corresponding Fokker-Planck equation (FPE) which has a nonlinear drift term 
and diffusion term with a time- dependent diffusion coefficient D(t). Our MD simulation reveals 
that the driven translocation process follows a SMperdiffusive law with a running diffusion coefficient 
D(t) oc t 1 where 7 < 1. This finding is then used in the numerical solution of the FPE which yields 
an important result: for comparatively small driving forces fluctuations facilitate the translocation 
dynamics. As a consequence, the exponent a which describes the scaling of the mean translocation 
time (r) with the length N of the polymer, (r) oc N a is found to diminish. Thus, taking thermal 
fluctuations into account, one can explain the systematic discrepancy between theoretically pre- 
dicted duration of a driven translocation process, considered usually as a deterministic event, and 
measurements in computer simulations. 
35 . In the non-driven case, / = 0, the translocation is slightly subdiffusive and can be treated within 

the framework of fractional Brownian motion (fBm). 

PACS numbers: 05.50.+q, 68.43.Mn, 64.60.Ak, 82.70.-y 
Ch ! I. INTRODUCTION 

o 
o 

The polymer translocation through narrow pores is one of the most striking single-molecule experiments motivated 
mainly by the DNA transport across cell membranes and the possibility of DNA sequencing [l|, In nearly two 
decades of intensive investigations considerable headway has been achieved by a variety of experimental, theoretical 
and simulational studies [3]. In most cases of experimental and simulation studies one deals with a polymer chain that 
^f-^ \ is driven by external electrical field which acts on the monomers passing through the nanopore from the cis- to the 
■ trans-side of the separating membrane [l[ . Some of the very basic questions one tries to answer concern the time r it 
I takes the chain to thread and translocate between the two half-spaces of the setup and in particular, how r depends 
I ■ on chain length N and applied force /. 
T — I ' The simplest theoretical approach has considered this process as one-dimensional biased diffusion (in terms of the 
fN) ■ translocated number of segments s) over an entropic barrier [J, Q ■ Such an approach obviously suffers from some 
inconsistency since the mean time of the driven translocation (r) turns to be smaller than the relaxation (Rouse) 
time of the chain tr so that the chain has not enough time to equilibrate and experience the global entropic barrier. 
Later, it was realized that the chain responds to pulling at first locally in the vicinity of the nanopore with a tensile 
force spreading away from the separating membrane. This scenario has been treated within the linear response theory 
(that is, for relatively weak driving forces) and the corresponding memory function has been derived explicitly @, 0] ■ 
For arbitrary strong driving forces, an interesting approach based on the notion of tensile force propagation along the 
chain backbone has been su gges ted by Sakaue [a-Hoj. Sakaue's idea has been used and worked out very recently in 
other theoretical treatments [ill [l2| . 

Despite increased insight in the nature of the translocation process, however, current understanding of the transloca- 
tion dynamics is still far from satisfactory. So, among other things, there exists a non-negligible systematic discrepancy 
between the theoretical (analytical) predictions concerning the main scaling exponents that characterize translocation 
dynamics and the data provided by computer experiments, mainly Molecular Dynamics (MD) simulation studies Q. 

Recently [HJ we suggested a theoretical description based on the tensile (Pincus) blob picture of a pulled chain 
and the notion of a tensile force propagation, introduced by Sakaue [8r1To| . Assuming that the local driving force is 
matched by a drag force of equal magnitude, (i.e., in a quasi-static approximation), we derived an equation of motion 
for the tensile front position. This enables one to calculate the deterministic dependence (i.e., without taking into 
account fluctuations) of the translocation coordinate as a function of time M(t). One can obtain the scaling law for 
the mean translocation time (t) vs. chain length N, i.e., (r) oc N'\ where a is the translocation exponent. The 
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total translocation process consists of two stages. First, the tensile force propagates so that more and more polymer 
segments get involved in the moving domain. Depending on the driving force /, which is applied to the segment in the 
pore, the moving domain attains different shapes: "trumpet" , "stem-trumpet" , or a "stem" , for weak, intermediate 
and strong forces, respectively. Once the last polymer segment on the cis-side of the pore gets involved, the velocity 
of the moving domain approaches a stationary value. After that the second stationary stage sets in and the rest of 
the chain is sucked into the pore with constant velocity. 

We have shown [l2[ that for the "trumpet" case, the characteristic duration (in the case of Rouse dynamics) of the 
first and second stages goes as n oc N 1+u / f and T2 oc N 2u / f 1 ^ respectively (where v 3/5 is the Flory exponent), 

so that at fn = aN v f jkoT 3> 1 the first stage time dominates. As a consequence, with growing chain length N and 
force /, the translocation exponent a also increases from a = 2v w 1.19 to a = 1 + v 1.59. These predictions are 
supported b y ou r simulation findings [l2| as well as by results of Lehtola et al. [1343 but differ from the results of 
Luo et al. Ilal . The tensile force propagation phenomenon treated within the iso-flux trumpet model also leads to 

a = 1 + v [njf 

Unfortunately, the MD-simulation results yield systematically smaller values for the translocation exponent a. For 
example, in our MD-simulations a ~ 1.33 for strong forces and a s» 1.06 for weak forces (distinct from the theoretical 
predictions a = 1 + v » 1.59, and a = 2v w 1.18 respectively). It is important to determine the origin of this 
inconsistency since apparently there is something missing in the afore mentioned theoretical consideration which gives 
room for speculations. For example, in a paper by Ikonen et at. [Tlj ] . the model based on the idea of tensile force 
propagation and the role of pore-polymer friction has been numerically investigated. The authors argue that 

the theoretical value for the exponent, a = 1 + v, may be seen only for very long chains whereas for the chain lengths 
used in real experiments or simulations the effective exponent a could be approximately 20% smaller. 

In this work we consider theoretically and by means of MD-simulations how fluctuations of the translocation coordi- 
nate (such fluctuations have been ignored in the quasi-static approximation as in the other deterministic treatments) 
affect the forced translocation dynamics. Indeed, the role of thermal fluctuations is by no means self-evident. A re- 
cent publication [l8| argues that only fluctuations related to the initial distribution of segments are essential whereas 
thermal fluctuations have a minor effect. Moreover, the authors argue [18| that the translocation exponents are not 
affected by the fluctuations in the initial conditions. 

In the present investigations we demonstrate that thermal fluctuations may facilitate the translocation dynamics 
so that the effective translocation exponent a becomes smaller. To this end in Section [II] we derive a Fokker-Planck 
equation (FPE) for the translocation coordinate s probability distribution function W(s, t). This FPE, which contains 
a nonlinear drift term and a time-dependent diffusion coefficient D(t), is then solved numerically demonstrating a 
significant fluctuation-induced facilitation of the translocation dynamics. In Section |TTT] we then present our simulation 
results. First, we prove that the translocation velocity under differently strong driving forces is a random process 
with a Gaussian distribution which is the only conjecture we need for our theory. Then we calculate the velocity 
autocorrelation function (VAF) which reveals an oscillatory behavior and a long-time tail. This gives rise to a 
time-dependent diffusion coefficient D(t) oc t 1 , where the exponent 7 < 1, indicating a swper-diffusively driven 
translocation. We end this report with a brief summary of conclusions in Section HVl Some technical details relegated 
to the Appendices E] [Bj 



II. FROM LANGEVIN TO FOKKER-PLANCK EQUATION 

In order to allow for fluctuations of the translocation coordinate one should consider the problem within the 
framework of the corresponding Fokker-Planck (FPE) equation which governs the probability distribution function 
(PDF) W(s,i) that exactly s segments have passed the pore at time t. By so doing we follow our earlier approach 
ref. 19] where the case of undriven translocation was treated. We focus here on the general case when the drift term 
in FP-equation is essential 

Our approach is based on the assumption that the translocation coordinate s(t), which measures the number of 
segments heaving reached the trans-side at time t, is a random process governed by the so-called V-Langevin (where 
V stands for "velocity" ) equation 

^ =»(•(*)) (2-1) 

This equation has been discussed extensively by Balescu (20j mainly in the context of plasma dynamics. The translo- 
cation velocity v (s(t)) in Eq. (|2.1[) . which depends on the trajectory s(t), is assumed to be a Gaussian process with 
given first moment 

(v(s(t))) d ^ K(s)\ s=s{t) ^ -1 U'(s)\ s=s{t) (2.2) 
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and a two-point correlation function 

(K S (*i)) - (v(s(h)))} [v(s(t 2 )) - (v(s(t 2 )))}) d = f G(h,t 2 ) (2.3) 

In Eq. (|2.2j) the function U(s) stands for the effective potential in the translocation coordinate space s while £o is 
the corresponding friction coefficient. The form of the averaged velocity K(s) as well as the effective potential U (s) 
we will find below using the correspondence principle. This principle states that in the limit of very small fluctuations 
the FPE solution reproduces the deterministic (or quasi-static) solution (s(t)) = M(t) which was derived in ref. [Hj]. 
The case of colored Gaussian statistics as well as the behavior of the VAF G{t\, t 2 ) will be studied by MD-simulation 
in Sec. IIIII This will be considered as the required dynamic input while solving the corresponding FPE. 

Based on the V-Langevin Eq. (|2.1[) . one could derive the corresponding FPE for the PDF W(s,t). The latter is 
defined as follows 

W(s,t) = (S(s-s(t))) (2.4) 

where (...) means an average over the Gaussian velocity field fluctuations with a first moment and two-point corre- 
lation function given by Eq. (|2.2I) and Eq. (I2.3[) . respectively. The detailed derivation is given in Appendix [S] The 
resulting FPE reads 



dW(s t) d d 2 

—m = -g- s ms)W( S ,t)]+D(t) — 2 W(s,t) (2.5) 
where the time dependent diffusion coefficient 

D(t) = f G(t,r)dr. (2.6) 
Jo 

It should be stressed that Eq. (12.51) is rather general. Indeed, no assumptions regarding, for example, the validity 
of the fluctuation-dissipation theorem are made in the course of derivation in contrast to the case of the generalized 
Langevin equation approach |2ll.l22j|. Moreover, the drift and diffusion terms are totally independent (i.e. the Stokes- 
Einstein relation does not necessarily hold). In other words, this FPE is valid even for a highly non-equilibrium driven 
translocation processes. 



A. What is the form of K(s)7 

As mentioned above, the form of K(s) can be fixed by making use of the correspondence principle which states that 
in the limit of zero fluctuations the FPE reproduces the deterministic solution for the first moment, (s(t)) — M(t), 
derived in ref. [Hj], i.e. 

M(t) = (cofarf (2.7) 

where the dimensionless force f a = af /T and time t — t/r with a, T and t = a 2 £o/T standing for the Kuhn segment 
length, temperature, and the microscopic characteristic time, respectively. In Eq. (|2.7I) cq is a constant of the order 
of unity and the exponent (3 in the limit of strong driving forces reads /3 = 1/(1 + v) [12(. 
In this limit G(t, t) = 0, D{t) = and Eq. ([23]) reduces to 

On the other hand, in this case W(s,t) = 5(s — M(t)) and Eq. (|2.8p is equivalent to 

p(s-M(t)) = -K(M(t))^- s 5(s-M(t)) 

= -M^-5(s-M(t)) (2.9) 

where in the first line we have used that K(s)5(s — M(t)) = K(M(t))5(s — M(t)), and in the second line the chain 
rule has been used. As a result the FPE is equivalent to the deterministic equation 

jM{t) = K(M(t)) (2.10) 
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It is easy to show that in order to restore the solution, Eq. (|2.7[) . the function K(s) in Eq. (I2.10p should be of the 
form 



at,o 



(2.11) 



It is also instructive to reconstruct the form of the effective potential U(s), given by Eq. (|2.2[) . Namely, from Eq. 
(I2.2[) and Eq. (I2.11[) one has U'(s) — —(/3cof/a) s^C 1- ^)// 3 . As a result, the effective potential that governs the 
translocation reads 



tt( \ tt P c °f 2-1/ 

U{S) = U °- a(2-l/P) S 

i.e., it has the form of the slanting plane shown in Fig. [1] 



(2.12) 
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FIG. 1: Effective potential U(s) 

fs>h> fi. 



given by Eq.( I2.12|l 'l as a function of the translocation coordinate s at different forces: 



B. Two ways to calculate D(t) by using MD-data 

The behavior of the diffusion coefficient D(t) which is necessary as an input for the solution of the FPE, could 
be found using a MD-simulation by two methods. First, one could directly calculate the VAF G(t, t) and then 
integrate it over the time according to Eq. (|2.6I) . On the other hand, there is a simple relationship between the 

time dependent variance (As 2 (£)) = f {[s(t) — (s(t))] 2 ) and D{t). Indeed, first of all s(t) = sq + J* v(s(ti))dti , i.e. 

(s(i)) = s + fo{v{s(t 2 )))dh. Thus, 



<[«(*) - M*)>] 2 ) = dh dt 2 (6v(s(h))6v(s(t 2 ))) 
Jo Jo 

= f dh f dt 2 G(h,t 2 ) (2.13) 
Jo Jo 



where, as before, Sv(s(t)) = v(s(t)) - (v(s(t))}. Finally, differentiation of (As 2 {t)) d = ([s(t] - (s(t)}} 2 } in Eq. (|2~T3)l 
leads to 



d(As 2 (t)) 



2dt 



f G(t,r)dT = D(t) (2.14) 
Jo 



where we have used that G(ti,t 2 ) — G(t 2 ,ti). In Section IIIII we will show that the mean squared displacement 
^As 2 (£)), calculated by the double integration of the VAF (see Eq. (|2.13p ) over time, and by the direct simulation, 
gives closely matching results. 
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C. Numerical Solution of the FP-equation 

The FPE, Eq. (|2.5I) . can be easily solved numerically assuming reflection-adsorption boundary conditions and some 
reasonable guess about the form of time-dependent diffusion coefficient D(t). The direct inspection of simulation 
data, given in Sec. IIIIB1 shows that the VAF G(ti,t2) has a long-time tail, i.e., correlations persist in time. This, 
according to Eq. (|2.14[) . leads to a time-dependent diffusion coefficient D{t) which could be approximated by a power 
law 

D(t) = — f> (2.15) 

TO 

where do is a constant and the exponent 7 < 1. Before we proceed further with the results of numerical solution, 
let us give a short discussion of the small- noise expansion of the FPE, Eq. (|2.5[) . where the diffusion coefficient is a 
small time-dependent function, i.e., D(t) — (e 2 /2)D (t), with D (t) cx t 1 . Following the book of Gardiner [23j . in the 
Appendix [B] we give a more extensive discussion of the small-noise expansion around the deterministic solution. In 
particular, for the first moment this expansion yields 

(#=(^V^?^ (2.16) 

As is evident from Eq. (|2.16[) . the fluctuations are responsible for speeding of the translocation process. The corre- 
sponding fluctuation increment is of the e 2 -order (the same as the diffusion coefficient, D(t) = e 2 Do(t)/2, itself) and 
goes with time as i 1+7_ ^. In order to estimate the exponent 1 + 7 — (3, one needs the exponent 7. One could use the 
data for (As 2 (i)), given by Bhattacharya et al. [24[ where it was reported that (As 2 (i)) ~ t 1A4 . As a result, accord- 
ing to Eq. (|2.14[) . the (reduced) diffusion coefficient Do(t) ~ t 0A4 , i.e. 7 = 0.44. This estimate leads (taking also into 
account that /3 = 1/(1 + v) « 0.63) to the behavior of fluctuation increment ~ i 81 which facilitates the translocation 
dynamics as (s(t)) — (co/i) ' 63 + (e 2 rfo// ' 63 ) i 81 . In the general case, an effective exponent enhancement is possible, 
if at least 1 + 7 — /3 > j3, i.e., 

7 > 2/3-1 (2.17) 

Thus, there is a lower limit for the 7 value which ensures the effective exponent (in the (s(t)) vs. t dependence) 
enhancement. For example, if (3 = l/{l + v) then 7 > (1 — v)/(l + v) « 0.26. 

The relationship, Eq. (|2.17l) . can be also readily obtained by the following arguments. If (As 2 ) oc t , then the 
first two moments could be estimated as (s) oc t e / 2 and (s 2 ) cx t e . But in the case of quasi-static (no fluctuations) 
approximation (s) cx t& and, provided that 8/2 > j3, some slope enhancement owing to fluctuations takes place. 
Taking into account that according to Eq. (|2.14l ) = 7 + 1, one immediately arrives at Eq. (|2.17l) . 

Figure [5] demonstrates the numerical solution of the FPE, Eq. (|2.5p . for the first statistical moment, (s(t)), as a 
function of time (the chain length N = 100). It can be seen that for a relatively strong force (upper set of curves for 
/ = 10), the fluctuation has little or no effect on the (s(t)) behavior. For weaker force, e.g., for / = 1, shown by the 
lower set of curves, and relatively large exponent 7 = 0.8, one can see a clear slope enhancement, up to 0.75. Thus, 
one may argue that under relatively small driving forces fluctuations facilitate the translocation dynamics. 

The results of numerical solution for the PDF W(s,t) as well as the corresponding statistical moments, (s(t)) , 
(s 2 (t)), (As 2 (i)), are shown in Fig. [3] for different parameters. As one might expect, the time dependence of the 
variance follows the law (As 2 (t)} oc t e , where the exponent 8 = 7 + 1 in accordance with Eq. (|2.14[) . In other 
words, for the running diffusion coefficient, D(t) oc i 7 , with 7 < 1, the exponent 1 < 8 < 2 and one finds a case 
of superdiffusion. It is of interest that, due to the external force (in this case / = 1) and the adsorption boundary 
condition (W(s — N,t) = 0), the PDF W(s,t) becomes narrower at a later stage of the translocation. This results in 
a non-monotonic behavior of the variance (As 2 (t)) for 7 = 0.4, and especially for 7 = (see Fig. [3]). These numerical 
findings for the statistical moments are qualitatively consistent with our MD-results given in Fig. [S] Eventually, the 
first passage time probability distribution (FPTD), Q(t) = —(d/dt) J Q W(s,t)ds, for different forces (/ = 1,5, 10) is 
depicted in Fig. SJ 

D. Continuous Time Random Walks (CTRW) approach and beyond 

It is pertinent to note that the interpretation of the superdiffusion in terms of the time persistent correlations in 
the VAF and the time-dependent diffusion coefficient goes actually beyond the polymer translocation problem. Recall 
that there are two reasons for the occurrence of an anomalous diffusion (25| : i) The consecutive steps of a random 
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FIG. 2: The result of numerical solution of Eq. (|2.5I) for two different driving forces, / = 1 (lower curves), / = 10 (upper curves), 
and exponents 7 = 0, 0.4, 0.8. The dashed line corresponds to the quasi-static approximation with a slope /3 — 1/(1 + v) ~ 0.63. 
Taking into account fluctuations with D(t) oc t°' s leads to a slope (or, an effective exponent) enhancement up to 0.75. The 
chain length N = 100. 



0.05 



0.02 



0.01 



(a) 
t=10 



t=ioo io 




10 10 10 t 10 
t=500 



10 

1=1000 



0.05 



0.02 ■ 



0.01 



20 



80 




1=2000 



FIG. 3: (a) The PDF function W(s,t) for / = 1, 7 = 0.4 and (inset) the corresponding statistical moments (s(t)), (s 2 (t)), 
(As 2 (t)} . (b) The same but for the exponent 7 = 0. 
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FIG. 4: The first-passage time distributions (FPTD) Q(t) calculated for 7 = 0.4 and different forces / = 1, 5, 10. 
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walker are independent but the waiting time distribution is a sufficiently broad function, i/j(t) ~ l/t 1+M , where 
fi < 1, so that the first statistical moment does no< exist ii) There is a long-time correlation between random steps. 
In the first case, the diffusion process could be rationalized in terms of a continuous-time random walks (CTRW) 
approach which finally leads to a subdiffusive behavior [2(1 1271. The anomalous diffusion in an external force could 
be quantified in terms of the fractional FP-equation (FFPeY|26|. [27| . Later, this FFPE formalism has been applied 
to the driven polymer translocation problem with an external potential approximated by a function linear in s, that 
is, U(s) = Uo — fs 28, 29]. Apparently, this approach does not take into account the tensile force propagation and 
the moving domain reflecting the chain reaction which were discussed above (see e.g. the resulting effective nonlinear 
potential Eq. (|2.12l) ). The calculations within the FFPE formalism suggest that to a leading order (s) cx and 
(s 2 ) oc t 2f * i.e. the variance goes superdiffusively, (As 2 ) oc £ 2m , because 1 < 2fi < 2 [29j . 

Our simulation results make it clear (see Sec. lIII Bl and Fig. (7^,) that in the case of translocation dynamics correlations 
persist in time and the CTRW-approach presumably could not be used. In contrast, the formalism of the V-Langevin 
equation, used in this Section, gives a general way to treat the anomalous diffusion. This approach corresponds 
closely to driven dynamics because neither the fluctuation-dissipation theorem or the Stokes-Einstein relation, nor 
time-translation invariance of the correlation functions appears to hold in this case. 



III. SIMULATION RESULTS 



A. Model 



In order to verify that the assumptions made in the theoretical analysis are justified, we performed a number of 
simulations of polymer chains threading through a nanopore. We used MD simulations as in Ref. fl9l ] . Here we briefly 
recapitulate the used algorithm. 

The model we used describes Langevin dynamics of a polymer chain which consists of N beads that thread through 
an octagonal pore in a closely-packed wall (membrane) . The interaction between the monomers of the chain is modeled 
by a Finitely Extensible Nonlinear Elastic (FENE) springs corresponding to a pair potential 

U F ene{u 3 ) = In ^1 - , (3.1) 

where r^ is the bond length between two beads and Rq — 1.5 is the maximal bond length. All beads experience 
excluded volume interactions which arc modeled by the repulsive part of the shifted Lcnnard-Jones potential, also 
known as the Weeks-Chandler-Andersen (WCA) potential. This potential Uwca is defined by 



UwcA{Uj) = 4e 



/ \ 12 








\nj] 





9(r c -r), (3.2) 



where O(x) is the Heaviside-function, i.e., we use a cut-off r c — 2~ 1 / 6 er, implying Uwca — for > r c . The 
monomers residing inside the pore experience a constant external force / in the direction perpendicular to the mem- 
brane, which we designate by x. The external force can be implemented by adding a linear potential U ex t, whose 
value is outside the pore and fx, if x is inside the pore region. Thus, / pulls the chain towards the region of positive 
x which we refer to as the irons-side. The equation of motion for the beads of the chain reads 

d?V' dv ' 

m—± = -V (Ufene + U LJ + Uext) - 7— 1 + Ki(t), (3.3) 
dt z dt 

where TZi(t) stands for a Brownian random force whose moments obey (TZi(t)) = and (TZf (ti)lZj (ta)) = 2fcsT7<5(ii — 
t2)fia,@$i,j The parameter values were set to e = 1.0, a — 1.0, k = 60.0, 7 = 0.70 and were kept fixed during the 
simulations. The temperature had a constant value given by T = 1.20e/fcs. 

The membrane is modeled by a plane of beads whose positions are kept fixed. Eight neighboring monomers are 
removed to obtain an octagonal pore. The interaction between the beads of the chain and the plane is mediated 
through the repulsive WCA potential. All simulations were performed starting from a configuration in which initially 
the chain is placed such that all but one monomer reside on the cis-side. In order to prevent the chain from escaping 
to the trans-side, we impose reflecting boundary conditions on the first monomer. The chain is fully translocated 
when all beads have made their way to the trans side. 
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FIG. 5: The velocity distribution, P(v z ), for two different forces, (a) / = 0, (b) / = 1, and three different time moments. Solid 
lines represent the Gaussian distribution fitting. 

B. Results 

First of all we need to prove that the translocation velocity follows indeed a Gaussian distribution. Simulations 
were performed for a chain length TV = 100 in the case of free (/ = 0) and driven (/ = 1) translocation. We have 
approximated the translocation velocity, v(t) — ds{t)/dt, by v z (t) , the z-component of the Cartesian velocity of the 
bead inside the pore [l9j . Then we recorded data of the velocity at specific time moments and made histograms over 
4000 runs. As can be seen from Fig. [5j the velocity corresponds closely to a Gaussian distribution (solid lines). This 
justifies our basic assumption, used in the derivation of the FPE, see Section IIT1 

The statistical moments of the translocation coordinate which are plotted in Fig. [5] for different forces, also 
provide an important information. While in the unbiased translocation case, Fig. [6^, the diffusion is Brownian (or 
slightly subdiffusive [Hj]), in the biased regime the process becomes superdiffusive, i.e., the variance (As 2 ) oc t e where 
1 < 9 < 2. Moreover, the exponent 9 increases with the growth of the driving force /, namely, 9 = 1.5 for / = 1 
(see Fig. |6Jd), and 9 — 1.84 for / = 10 (see Fig. |6jl). For relatively large forces (/ = 5, 10), the variance (As 2 ) is 
nonmonotonic, i.e., it goes through a maximum at a late stage of the translocation. This finding is in agreement with 
our results of the FPE numerical solution given in Sec. Ill CI (see Fig. [3|). 

The superdiffusive behavior for (As 2 ) implies (according to the relationship Eq. (12.14)) ) that the diffusion coefficient 
is a growing function of time which has been approximated by a power law, Eq. (|2.15[) . in Section [II CI On the other 
hand, in accordance with Eq. (|2.6j) . the time-dependent diffusion coefficient D(t) can be obtained by integration of 
the VAF G(t 1 ,t 2 ). In Fig. [TJi we display the normalized VAF G(ti,t 2 )/G(t 1 ,t 1 ) for t 1 = as a function of t 2 . This 
figure indicates that the VAF has a fairly complex behavior: it changes sign and reveals a long-time tail. It can also 
be seen that stronger forces imply longer time correlations which give rise to the observed superdiffusive behavior. 
All measurements have been averaged over 5000 runs (except for / = 1 which was obtained by averaging over only 
400 runs). 

The corresponding diffusion coefficient D(t) could be obtained by simple integration over time, according to Eq. 
(I2.6|) . The result of this calculation, given in Fig. [TJd, shows that the exponent 7 attains different values within 
different time intervals. Eventually, we perform a consistency check. On the one hand, we integrate the expression for 
G(ti,t 2 ) over time arguments which provides, according to Eq. (|2. 13|) . the value of (As 2 (t)) = J Q dt\ J dt 2 G(ti 7 t 2 ). 
On the other hand, the MD-simulation yields the direct time dependence (As 2 (i)). This comparison of two different 
calculations of (As 2 (t)), displayed in Fig. [SJ gives almost identical results for / = 10. 

IV. CONCLUSION 

We have investigated theoretically as well as by extensive MD-simulation the role of fluctuations for the case of driven 
polymer translocation through a nanopore. The consideration is based on the V-Langevin equation, Eq. (|2.1[) where 
the velocity is a random process with Gaussian distribution. Keeping this in mind we have derived a corresponding 
FPE, Eq. (|2.5p . with a nonlinear drift-term and time-dependent diffusion coefficient D(t) which can be represented as 
time integral of VAF G(ti,t 2 ) (see Eq. (|2.6p ). The derivation requires neither the validity of the fluctuation-dissipation 
theorem (as in case of the generalized Langevin equation approach [2l| ) nor the time-translation invariance of G(t 1 , t 2 ) . 
Such a general approach becomes an important tool for the description of driven translocation dynamics. 

It is pertinent to note two limiting cases. In the zero-fluctuation limit, e.g., when the driving force is large, one 
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FIG. 6: The first (s) and second (s 2 ) moments as well as the variance (As 2 (t)} — (s 2 ) — (s) 2 for chain length N — 100 and 
different driving forces /: (a) / = 0. (b) / = 1.0, (c) / = 5.0. (d) / = 10.0. 
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FIG. 7: (a) The velocity autocorrelation function for four different values of force. The stronger forces trigger a longer in time 
correlations. The vertical axis is normalized by G(ti,ti). (b) The diffusion coefficient D(t) as function of time. Larger diffusion 
coefficients arise with / increases. 



recovers the deterministic solution. For example, for the first moment it is given by Eq. (12.71) ) which we have discussed 
recently [l2j . In the zero- force case (nondriven translocation), in contrast, Eq. (|2.5p leads to the fractional Brownian 
motion (fBm) description which was treated in our previous paper fl9j|. The fact that D(t) is a growing function of 
time (see Eq. (|2.15t "i . implies that VAF G(ti,t2) has a long-time tail (long persistence of correlations). This has been 
verified by our MD-simulation results (see Fig. [7^, and Fig. [7}j ). By making use the MD-simulation we have also 
checked the Gaussian distribution of the translocation velocity (see Fig. [5j . 
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FIG. 8: The value of (s 2 (t)) calculated by two different ways: by integrating the VAF G(ti,t2) over time, i.e. (As 2 (t)) = 
fg dt\ f Q dt2G(ti,t2), and by direct simulation. Chain length N = 100, force / = 10. 



The numerical solution of the resulting FPE, Eq. (|2.5I) . reveals a number of salient features. If the exponent 7 
governing the time dependence of the diffusion coefficient, D(t) oc V , is large enough (see the condition given by 
Eq. (|2.17j) ), fluctuations "assist" the translocation dynamics and the translocation coordinate (s) oc dependence 
becomes steeper. Moreover, the variance follows the law (As 2 ) oc t e , where the exponent 1 < 9 < 2, i.e., the dynamics 
is superdiffusive. The increase in slope due to fluctuations could, therefore, explain the systematic disagreement 
between deterministic theory and the MD-simulation results which has been discussed in ref. jl2j . One should recall 
that the scaling law for the mean translocation time reads (r) oc N a where a w 1//3. Owing to fluctuations, (3 grows 
so that the translocation exponent a becomes smaller in agreement with the MD-simulation. 
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Appendix A: Derivation of the Fokker-Planck-Equation (FPE) 

Differentiation of Eq. (12.41) with respect to t and using the chain rule and the property of the (^-function yields 

dW(s,t) 8 /ds(t) 

-S(s-s(t)) 



dt ds \ dt 

= ~ <«(*(*)) *(«-«(*))) (Al) 

where the differential operator d/ds can be put out of averaging because the S- function is the only one which depends 
on s. In Eq. (|A1| we have used the Langevin equation Eq. (|2.1|) . By making use the relation (v(s(t)))5(s — s(t)) = 
{v(s))5(s — s(t)), Eq. (jAip can be written as 

= ~K(s) (8(a - *(*))> - (6v( 8 (t))8(8 - s(t))) (A2) 

where dv(s(t)) = f v(s(t)) — (v(s(t))) and we have used Eq. (j2.2j) . The second term in Eq. (|A2I) could be expressed 
in terms of W{s,t) by employing Novikov's theorem [30L l3l|. According to Novikov's theorem, if g{t) is a colored 
Gaussian random process with a zero average, i.e., (g(t)) = 0, and the two-point correlation function is given by 
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(.9(^1)5(^2)} = G(t\,t2), then for an arbitrary functional, lZ{g(t)}, the average (g(t)lZ{g(t)}) can be written as 

t 

( 9 m{g(t)}) = J drG(t,T)( Sn ^ } ^ (A3) 


where the symbol S/Sg(r) stands for a functional derivative. 

Let's use now the Novikov's theorem to recast the second term in the r.h.s of Eq. (|A2[) . In this case g(t) = 8v(s(t)) 
and lZ{g(t)} = 8(s — s(t)). As a result, Eq. (IA2[) takes on the form 

9W(s,t) _ 3 r __ d } 8 s(s _ s{t)) 



.tj-[K(s)W(8,t)]-^ J dr G(t,r) 



dt ds ds J \5v(s(t)) 


t 

= -°. m e) W{ s,t)] + * i J *G{t,r)(6(s-sit))^) (A4) 



where the random variable s(t) is treated as a functional of v(s(t)). But s(t) = s(0) + f Q v(s(t'))dt' and one has 



because r <t. As a result, one obtains the FPE Eq. (|2.5p ~ (|2.6[) . 

Appendix B: Small noise expansion 

In this Appendix, following Section 6.3 of the book of Gardiner (23|, we give a short exposition of the small- noise 
expansion for the FPE, Eq. (|2.5p where the time-dependent diffusion coefficient D(t) = e 2 Do(t)/2, with Do(t) = 

(d /T )t-y. 

dWis t) d e 2 d 2 

= _£_ [K{s)w(Sit)] + z_ Do{t) WM ( Bi) 

In the same way as in [23| . one can expand around the deterministic solution, i.e., 

s = M{t)+ex (B2) 

where a; is a new random variable. Here we discuss only the first statistical moment which is measured by our 
computer simulation experiment and has the following form (to the order e 2 ) (see Eq. (6.3.16) in |23j) 

<«(*)> = M(t) + e (x{t)) = M(t) + eXl{t) + e 2 X\{t) (B3) 

where the evolution of X${t) and X\(t) is given by three ordinary differential equations (see the corresponding 
Eqs. (6.3.19), (6.3.20), (6.3.22) in (H) i.e., 



dt 
dX\{t) 

dt 
dX 2 (t) 

dt 



K 2 (t) X 2 (t) + K x {t) X\{t) 

2K t (t) X 2 (t)+D Q {t) (B4) 



where K\(t) and K^it) are two expansion coefficients of K(s) around the deterministic solution, i.e., K\(t) 
\dK{s) / ds] s -M{t) and Ki(t) — (l/2)[d 2 K(s)/ds 2 ] s= M(t)- Taking into account Eq. (|2.1ip . one arrives at 



T t 
1 (l-P 

to I 2(3^ J P 



*»(*) = - ( TTT I TflTTTs (B5) 
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where the dimensionless force and time are given as / = af /T and t — t/ro, respectively. 

In order to solve the first order, linear, inhomogeneous ordinary differential equations, Eqs. ( IB4|) . one recalls that 
the corresponding generic equation has the following form 



dt 



y(t) + a(t) y(t) = b(t) 



It easy to verify that this equation has a solution (see [32| Sec. 4.3) 



y(t) = exp{- / a(t)dt} 



exp{ / a(t)dt} b(t)dt + c 2 



(B6) 



(B7) 



where c 2 is a constant of integration. 

Now let's go back to Eqs. (|B4[) . Solution of the first equation in Eq. (|B4[) is given straightforwardly as 



Xl{t) = Cl exp{/ K^dt'} = ci exp{-(l-/3)ln(Vr )} 

J T 

= ^ (B8) 

where c\ is a constant. Then, taken into account Eq. (|B7[) and Eq. (|B5[) . the solution for the third equations in 
Eq. (|B4|) is given by 



r \ 2 ( 1 -' 3 ) 



t / t ^2(l- 



\ r o 



D (t') dt' 



(B9) 



It is pertinent to note that while using Eq. (|B7|) we set c 2 = because the fluctuation corrections cancel as soon as 
the diffusion coefficient is zero, i.e., X§(t) = at do = 0. 

Again, by making use Eq. (|B7[) as well as of Eqs. (|B5|) and (IB9|) . the solution of the second equation in Eq. (|B4|) 
reads 



(1-/3) 



4 'i'^ 1 " 



\ T o 



K 2 {t') Xl{t')dt' 



JO /TO 



//3 V t 



(1-/3) 



V r 



(1-/3) 



(1+/?) Cl+7) 
TO 



dt 1 

TO 



(BIO) 



Collecting all results, given by Eqs. (|B3]l , (|B8]l and f(BTo|) . one has 

«'» = (^)^A + £ F 

which goes to Eq. (|2.16p for large times, i.e., t > 1. 



ii+7- 



(Bll) 
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